Radio quiet neutron star IE 1207.4-5209 as a source of gravitational waves 
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There are four puzzles on IE 1207.4-5209: (1) the characteristic age of the pulsar is much higher 
than the estimated age of the supernova remnant; (2) the magnetic field inferred from spin-down is 
significantly different from the value obtained from the cyclotron absorption lines; (3) the spinning 
down of the pulsar is non-monotonic; (4) the magnitude of the frequency's first derivative varies 
significantly and its sign is also variable. The third puzzle can be explained by a wide binary system, 
with orbital period from 0.2 to 6yr. This letter proposes that all four puzzles can be explained 
naturally by an ultra-compact binary with orbital period of between 0.5 and 3.3 min. With the 
shortest orbital period and a close distance of 2kpc, the characteristic amplitude of gravitational 
waves is h ~ 3 x 10~ 21 . It would be an excellent source of gravitational-wave detectors such as the 
Laser Interferometer Space Antenna. 

PACS numbers: 04.30.Db, 97.60. Jd,98.70.Rz 



I. INTRODUCTION 

The radio quiet neutron star IE 1207.4-5209 (here- 
after 1E1207) is at the center of the supernova remnant 
(SNR) PKS 1209-51/52. It was discovered by Helfand & 
Becker 1] with the Einstein Observatory. The distance 
to the SNR is d = 2.1 +1 Q % kpc @. The X-ray spectrum of 
the central source can be described by a thermal model 
which gives a distance of 2 kpc 0, |j, • 

The long observations devoted to 1E1207 both by 
Chandra and by XMM-Newton have unveiled a num- 
ber of unique and somewhat contradictory characteristics 
that, at the moment, defy standard theoretical interpre- 
tations [||. 

The characteristic age of the pulsar, 200 to 900 kyr 0, 
is much larger than the estimated age of the SNR, 3 to 
20kyr Q. 

The values of the spin-down luminosity, E ~ 1 x 
10 34 ergs _1 , and conventional magnetic field (B-ficld), 
£?~3xl0 12 G, are typical for a radio pulsars How- 
ever such a B-ficld is significantly different from the value 
obtained from the cyclotron absorption lines interpreted 
both in terms of electrons (B ~ 8 x 10 10 G) as well as 
protons (B ~ 1.6 x 10 14 G) [13, EU 

Chandra and XMM-Newton observations indicate that 
the pulsar is not spinning down steadily 0,0,0] . More- 
over, the first derivative pulse frequency varies signifi- 
cantly and its sign is also variable in different observa- 
tions 

The non-monotonic behavior of its pulse frequency, 
Az/, is interpreted by three hypotheses: glitch, accretion 
and binary (with an orbital period of 0.2 to 6 yr), in which 
the binary hypothesis is somewhat more plausible than 
the other two . 

This letter analyzes that for binary pulsars with very 
small mass function, the Roemer time delay in one or- 
bital period cannot be resolved. Such pulsars may thus 
be treated as "isolated" pulsars. However the pulse fre- 
quency and frequency derivatives of such pulsars are still 
affected by the orbital motion at long time scale, which 



causes anomalies like that of 1E1207. 

Different from the orbital period of 0.2 to 6yr 9], the 
one that predicted by this letter is much shorter, 0.5 to 
3.3 min, therefore, 1E1207 is an ideal source of gravita- 
tional waves. 



II. ORBITAL EFFECT AT LONG TIME SCALE 

Roemer time delay is the propagation time across the 
binary orbit, which is given, 



r sin z 



sin(w + /) . 



(1) 



where c is the speed of light, r the distance between the 
focus and the pulsar, / the true anomaly, uj the angular 
distance of the periastron from the node, and i the orbit's 
inclination. The orbital motion also causes the change of 
pulse frequency, Ab>, 



K [cos(w + /) + ecosw] , (2) 



where K = 2ira p sini/[cPb(l — e 2 ) 1 ' 2 ) is the semi- 
amplitude, e, Pb, Op are eccentricity, orbital period, and 
pulsar semi-major axis, respectively. 

Small companion mass, i or Pb of a binary pulsar may 
make the Roemer time delay of Eq. Q unmeasurable. A 
binary pulsar may thus be treated as an "isolated" pulsar. 
Whereas, following calculation indicates that Eq. Q and 
Eq. J2J) can still cause long-term effects on such "isolated" 
pulsars. 

For a binary pulsar, the time received by the observer 
(barycentric time) is, 



tb = t P H — 



(3) 



where t p is the proper time of the pulsar, and z/c is 
dependent on Kepler equation, 



E — e sin E = M = nt . 



(4) 
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where M , E and n are mean anomaly, eccentric anomaly 
and mean angular velocity, respectively. Note that t is 
the time of periastron passage, which is uniform. 

For a true isolated pulsar, we have z/c — in Eq. (J2J), 
thus ib = t pi which means both t\, and t p are uniform. 
But for a binary pulsar system, tb is no longer uniform, 
whereas t p is still uniform. 

Therefore, the proper time of the pulsar, t p , can be 
used to replace the uniform time, t of Eq. (J3J, then we 
have M = nt p . 

If Av of Eq. is averaged over one orbit period by 
the measured time, t^, then it gives 



(At/) = — 



Audtu = — 
P 



b Jo 



Av(dt p + -dt p ) 

c 



1 f b z X f b 

— / Av-dt p = — / Av cos(w + E)Edt p 

Pb Jo c Pb Jo 

XKv 



n n(l--) + 0(e 4 ). 



(5) 



where X is the projected semi-major axis, X = a p sin z/c. 

In practical observation, an observer averages Av from 
to T (T ~> Pb) through t D , the time received by ob- 
server, without knowing the orbital period, Pb, at all. 
However if the pulsar measured is truly in a binary sys- 
tem, P b will affect the averaged result, as given in Eq. J5J, 
thus the averaged Av given by the observer is 



(M = - 



Avdth = — 



P b N 



Avdt h + 



Pb(N-l) 



Avdt 



P b N 




Avdty 



+ o((3f) , (6) 



where (3 = XKvtt(1 — e 2 /4)/P\>, and N is an integer. 
Eq. © indicates that if a pulsar is in a binary system, 
then (Av) measured by the observer is actually contam- 
inated by the long-term orbital effect, (3. 



III. INTERPRETATION OF FOUR PUZZLES 
AND ESTIMATION OF ORBITAL PERIOD OF 
IE 1207 

A. Puzzle 1: pulsar age vs SNR age 

If the Av of Eq. JSJ and Eq. © (brackets (, ) are ig- 
nored hereafter) are unchangeable then the effect actu- 
ally cannot be measured. However [3 contains orbital el- 
ements, i, e and a (where a is the semi-major axis of the 
orbit, a = a p M/M2, M and Mi are the total mass and 
companion mass respectively, the mass of the pulsar is 
Mi) which are long-periodic terms when the Spin-Orbit 
coupling effect is considered. Therefore, Av is a function 
of time, and the orbital effect induced !>l can be found 



by differentiating (3 of Eq. ©, as given in detail in the 
following subsection. 

Thus the observational first derivative of the pulse fre- 
quency, !> obs , is given by 



I'obs 



(7) 



where v is the intrinsic one, which is caused by magnetic 
dipole radiation. Thus the following relation can be ob- 
tained 



2t>obs _ ,2*> + 2v L v ^ 2;> ^ 2vl_ 
v ohs v v v ohs v v 



Eq. JHJ is actually 



2v L 



(8) 



(9) 



where r = 200 — 900 kyr is the age corresponding to the 
contaminated spin-down (by the long-term orbital effect). 
In other words, when the true age of the pulsar equals 
the age of SNR, r p = 3 — 30 kyr is the true characteristic 
age of the pulsar. 

Putting t and r p into Eq. one obtains two group 
solutions corresponding to maximum and minimum mag- 
nitude of vl and v, respectively, (vl = 1.2x 10~ n Hz s _1 , 
v = -1.3 x 10~ n Hz s" 1 ); and {v L = 1.7 x 10~ 12 Hz s" 1 , 
v = -1.9 x 10~ 12 Hz s" 1 ). 

This implies that the magnitude of v and vl are much 
larger than that of ^obs, since vl and v nearly cancel each 
other out. Therefore the age puzzle, r p <^ r, can also be 
explained. 



B. Puzzle 2: B-field 



Section flll Al shows that the measured !> b s may under- 
estimate the true intrinsic pulse frequency derivative, v. 
Therefore, the B-field 3 x 10 12 G inferred from i> b s may 
be under- estimated also. 

The two !>l obtained through Eq. correspond to 
two magnetic dipole radiation-induced v, and therefore, 
to two B-field, 3 x 10 13 G and 1 x 10 13 G, respectively. It 
is easy for them to reconcile with the high B-field option, 
B = 1.6 x 10 14 G 10], i.e., by assuming the magnetic 
inclination angle, a = 11° and a = 4°, respectively. 

However, it is very difficult for these two B-fields to rec- 
oncile with another option, B — 8 x 10 10 G 01- There- 
fore, the B-fields inferred from the true intrinsic spin- 
down favors that 1E1207 is a magnetar. 



C. Puzzle 3: non-monotonic spin-down and 
estimation of orbital period 

In the gravitational two-body problem with spin, each 
body precesses in the gravitational field of its compan- 
ion (geodetic precession), with precession velocity of 1 
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Post-Newtonian order (PN) |l7|. The Spin-Orbit cou- 
pling causes long-periodic variations in the six orbital 
elements, i, e, a, M, uj and (longitude of the ascending 
node) By the definition of K and X, (3 of Eq. © 

can be rewritten as, 



(10) 



where p = 7rsin i (M 2 /M) 2 (1 - e 2 /4)/\/l - e 2 . Accord- 
ing to Eq. Q10[l . /? contains the orbital elements, e, i and 
a, which are all long-periodic terms when the Spin-Orbit 
effect is considered. However, the variation of i is much 
smaller (S/L times, S and L are the spin and orbital an- 
gular momentum respectively) than that of a and e |l8j |. 
Thus the long-period variation of Eq. © can be written 
in a Taylor series as 



Av = = fa +fit + ... = fa ~ - (1 - 0* + - . (11) 

a 



where 



(1 - e 2 )e 2 



2(l + e 2 )(l-e 2 /4) 1 + e 2 



(12) 



and 



GP(l+e 2 ) 

c 2 a 3( 1 _ e 2)5/2 



(2- 



3M 2 
2Mi 



) (PyQz - PxQy) (13) 



where P x , P ; , C^, Q y are sine and cosine functions of 
w and f2 [lj . The orbital period Pb of a few minutes 
corresponds to Cogr ~ 10~ 5 s _1 , which corresponds to 
a/a ^ 10~ 6 s _1 . Thus in the observation time span, At ~ 
10 2 ks, a has changed like, one tenth of its period, which 
means u> has changed by 7r/5. This actually corresponds 
to a large variation amplitude in Ao. 

Define da = |Aa/a| max , then the maximum and min- 
imum a of Eq. H1(J|) are a max = a(l + <5a), and a m i n = 
a(l — <5a) respectively. The discrepancy in Av corre- 
sponds to the error bar of each observation is given, 



5v = 



1 



1 



GMvp 
2nc 2 a \ \ — 5a 1 + 5a 



25a 



1 - (5a) 2 



(14) 



The fact that the amplitude of 5v is not much larger 
than a few //Hz || demands that a > \Aa\. Thus the 
maximum 5a can only be like 5a = 0.9, whereas 5a — 1 
is not allowed. 

From the point of view of Eq. QlOjl. both error bars 
in one observation and discrepancy for different observa- 
tions are dependent of the variation of a. The difference 
is that the discrepancy among different observations, i.e., 
Jan 2000, Aug 2002, corresponds to a much longer time 
scale, in which Av is modulated by both u) and f2 (the 
period of Q is comparable to that of u>) for many pe- 
riods. Whereas in one observation (10 1 ks- 10 2 ks), the 
time may be just enough for Av to vary in a few periods 
of a;, or even less than a period of u>. 



The jump of Av between Dec 2001 and Jan 2002 0, 
can be explained by the variation of w and f2, which 
causes relatively sharp variation in a and thus significant 
variation in Av. 

The second term at the right hand side of Eq. 111|) 
actually corresponds to vl, which is given in magnitude 
as 



v L =p-(l-0 
a 



(15) 



Putting the two j>l obtained in Eq. @ into Eq. (|15fl . we 
have two curves p vs Pb corresponding to v^i and ii\jx 
respectively, as shown in Fig. ^ 

Similarly putting Av\ = 0.18 //Hz and Av 2 = 4.2 /iHz 
into Eq. J^J, we have two p vs Pb curves as shown in 
Fig.^ which correspond to minimum and maximum dis- 
crepancies in v or error bars in different observations of 
Zavlin et al ■ 

The maximum orbital period, 3.3 min, is given by the 
cross section of A v 2 and i>L2 at C as shown in Fig. ^ 
The minimum orbital period is 0.1 min corresponding to 
A given by Avi and i'Li- 

In the area ABCD of Fig. ^ a point, i.e., with 
Pb = 0.7 min and p = 0.008 can be found. Assuming 
Mi = 1.4M , M 2 = O.2M and e = 0, p w 0.09 sin 2 i 
is obtained, and by the definition of p, s'mi w 0.3 can 
be obtained. In turn X is given, X 2.6 ms, which is 
smaller than the time resolution of observation, 5.7 ms 
or 2.9 ms 0. Therefore, the modulation induced by the 
orbital motion may not be detected from the timing ob- 
servation. This is consistent with the fact that the side 
band corresponding to Pb of a few minutes has not been 
found in 1E1207. 

Therefore, the companion of 1E1207 should have low 
mass and be compact enough. Because the undetected 
orbital modulation implies that X must be small; and 
the short orbital period demands that the companion be 
compact object like low-mass neutron star 01 or strange 
star [TH Il6j |. A white dwarf star companion is unlikely 
due to the separation of the two stars is almost equal to 
the radius of a white dwarf star when the orbital period 
is of 1 min. 



D. Puzzle 4: magnitude and sign of z> b s 

The values of Pb, p and the period of uj corresponding 
to the four points ABCD of Fig.QJare shown in Tabled 

As given by Eq. (fTOfl and Eq. lfT3|) . both Av and i> L 
vary with a which is in turn modulated by the period, T w 
(and Tq) as shown in Table UJ In the case Pb = 1 min, 
the period of T w is ~ 25 days, therefore the period of 
variation of uj_, is approximately ~ 25 days also (recall 
the period, To, is comparable with T u ). 

On the other hand, the intrinsic v changes steadily, 
which means v and z/l some times cancelling out, and 
some times have the same sign and enhancing, thus v Q \, s 
can both be ~ 10 -14 s -2 (v and z>l cancelled out); and 
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TABLE I: Results correspond to four points of Fig. 



point Pb(min) 


P 




T,,, (day) 


A 0.1 


6 x 10" 


4 


0.6 


B 0.3 


1 x 10" 


3 


4 


C 3.3 


3 x 10" 


1 


183 


D 1.0 


1 x 10" 


1 


25 



To, — 2-Ka/a represents the period of a /a. 



~ 10~ n s~ 2 (z> and i> L enhanced). This well explains 
the observations of Zavlin et al which show that 
i'obs can have very different magnitude, ~ 10~ 14 s~ 2 and 
~ 10~ n s~ 2 and its sign is also changeable at different 
epochs. 

Eq. (|10[1 and Eq. I|15|> actually predict that Av and j>l 
can vary with periods of days, thus v hs can change sign 
in order of days, or even during one observation (10 2 ks). 
Comparing Table [I] with observation may extract the pe- 
riod T w and Tq and therefore determine the orbital pe- 
riod Pb- 

IV. DISCUSSION 

Therefore, all four puzzles can be explained naturally 
by an ultra-compact binary system. 

The best spectral model describes the continuum as the 
sum of two blackbody curves with kT = 0.211±0.001kev, 
for an emitting radius R = 2.95 ± 0.05 km; and kT = 
0.40 ± 0.02 kev (R = 250 ± 50 m) 13]. It is possible that 
these two emitting radii are from the hot spot of two 



stars, 1E1207 and its companion. 

The characteristic amplitude of gravitational waves 
from a binary system is [jjj 

~3xl0~ 21 , (16) 

where /i is the reduced mass which equals \x — 0.18A/ Q 
when Mi = 1.4M and M 2 = 0.2M©; d = 2kpc is the 
distance; Pb = lmin is orbital period; and /(e) is given 
by /(e) = (l + |fe 2 + §e 4 )/(l-e 2 ) 7 / 2 , which is assumed 
/(e) « 1. In such case, GM/(c 2 a) ~ (v/c) 2 - 9 x 10~ 5 
(v is a characteristic orbital velocity), which means equa- 
tions, Eq. Ullfl -Ea. iJTSJl, based on the Post-Newtonian 
approximation are good enough to describe the dynam- 
ics of the ultra-compact binary system. 

The time scale of coalescing of 1E1207 corresponding 
to Pb = 0.5 min is ~ 27 yr. In order to be consistent with 
the fact that it was discovered in 1984 and is still there, 
it is necessary that Pb > 0.5 min. Therefore the most 
probable orbital period of 1E1207 is (0.5-3.3) min. 

The low wave frequency, ~ 10 -2 Hz, and the extremely 
large wave amplitude means that 1E1207 is an ideal 
source for the space detector Laser Interferometer Space 
Antenna. 
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FIG. 1: The two Av curves (Aui = 0.18 ^Hz, and 
Av2 — 4.2 /iHz) intersect with two !>l curves (uli = 1.2 x 
10 _11 Hz s _1 and i> L2 = 1.7 x 1(T 12 Hz s _1 ) at four points 
(A, B, C and D). Therefore, the binary parameters of 1E1207 
are contained in the area ABCD, in which the minimum and 
maximum values of p (= it sin 2 i(M 2 /M) 2 (1 — e 2 /4) — e 2 ) 
are from 4x 10~ 4 to 2x 10 _1 , and those of Pb are from 0.1 min 
to 3.3 min. Considering the coalescing time constraint, the 
minimum Pb should be 0.5 min. The vertical arrow repre- 
sents Pb = 0.5 min. 



